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Quantum emitters coupled to surface plasmons of a nano-wire: A Green function approach 
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We investigate a system consisting of a single, as well as two emitters strongly coupled to surface plasmon 
modes of a nano-wire using a Green function approach. Explicit expressions are derived for the spontaneous 
decay rate into the plasmon modes and for the atom-plasmon coupling as well as a plasmon-mediated atom-atom 
coupling. Phenomena due to the presence of losses in the metal are discussed. In case of two atoms, we observe 
Dicke sub- and superradiance resulting from their plasmon-mediated interaction. Based on this phenomenon, 
we propose a scheme for a deterministic two-qubit quantum gate. We also discuss a possible realization of 
interesting many-body Hamiltonians, such as the spin-boson model, using strong emitter-plasmon coupling. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

A strong coupling between individual quantum emitters and 
photons is one of the key ingredients for photon-based quan- 
tum information processing. It is needed to achieve a reli- 
able transfer of excitation between stationary and flying qubits 
and to realize quantum logic operations. Strong radiative cou- 
pling requires a tight confinement of the electromagnetic field. 
This can be achieved e.g. when quantum emitters are coupled 
to surface plasmons of nano- wires as proposed in |Q]]. Be- 
sides its simple structure, the atom-nano-wire scheme guar- 
antees an exceptionally strong field-emitter coupling. The 
original setup involves a single atom placed close to the sur- 
face of a metallic wire having sub-wavelength radius. In the 
present paper we analyze this scheme in detail by means of a 
Green function approach taking into account metal losses and 
extend it to a plasmon-mediated coupling between different 
atoms. The unique properties of the propagating electromag- 
netic modes of the wire, the surface plasmons, have opened 
up new possibilities in many fields such as waveguiding be- 
low the diffraction limitpfl, sub-wavelength imaging y| or 
enhanced fluorescence |4[]. The small transverse-mode area 
also enables one to strongly couple an emitter to the plasmon, 
giving rise to a substantial Purcell effect with a Purcell factor 
of 10 2 or larger. The system has been proposed as an effi- 
cient single-photon generator |Q]] as well as a single-photon 
transistor ||5|], strong emitter-plasmon coupling has been ex- 
perimentally shown using quantum dots [6], single plasmons 
alongthe wire have been detected, using N-V centers as emit- 
ters QV|] and the wave-particleduality of surface plasmons has 
been verified in experiment Jit . 

We describe a system where a single, as well as two emitters 
interact with the surface plasmons of the wire in a fully quan- 
tum mechanical approach. Due to the strong coupling, the 
plasmons mediate an intensive long-range interaction between 
different emitters. As a result, we observe sub- and superradi- 
ance in a two-atom system which can be efficiently controlled 
by adjusting the distance between the atoms. Based on this 
phenomenon, we propose a possible scheme for a determin- 
istic quantum phase gate involving 3-level lambda atoms and 
classical 2ir laser pulses. 



The paper is organized as follows: In Section [TT] we intro- 
duce the Green tensor of the system as well as the structure 
and properties of the guided modes. Section [III] discusses 
the coupling of a single emitter to the plasmons. The total 
spontaneous decay rate and the decay rate into the plasmon 
modes are calculated. We also derive an expression for the 
atom-plasmon coupling strength. Section [IV] deals with the 
long-range interaction of two emitters coupled to the plasmon 
modes. We show the emergence of Dicke- sub- and superra- 
diance which can be controlled by adjusting the distance be- 
tween the emitters. We present a scheme for a deterministic 
quantum phase gate based on the controllable superradiance 
phenomenon. Both in the case of a single emitter, as well as 
a pair of them, we discuss the effects of metal losses and their 
consequences. Finally, we discuss a possible extension of the 
level scheme of the emitters to realize interesting many-body 
Hamiltonians such as spin-boson model using the coupling to 
nanowires. Section[Vlcontains a summary of the results. 



II. ELECTROMAGNETIC GREEN TENSOR AND GUIDED 
MODES OF THE WIRE 

The coupling of quantum emitters to the quantized elec- 
tromagentic field in the vicinity of linear macroscopic media 
can be described by incorporating the effect of the medium 
in the Green tensor of the field G(r,r',uj) which obeys the 
Maxwell-Helmholtz wave equation with a Dirac delta source 
at r = r ' 



V x V x -e(f, w) 

c 2 



G(r,f',u}) = IS(F-r'). (1) 



Here e(r, oS) is the complex dielectric function of the isotropic 
linear medium at position r and frequency uj. In the following, 
we will use e = — 75+0. 6i as typical value which corresponds 
to the electric permittivity of silver at around A = lpm. It is 
assumed that the magnetic response is negligible. 

In the following we will calculate the Green tensor for a 
nano-wire. We consider a cylindrical wire with radius R and 
symmetry axis z. The metal is described by a complex dielec- 
tric function with e' = Re[e] < 0. For simplicity we assume 
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an infinitly long wire to have translational invariance along the 
z-axis. For f' being outside the cylinder we can decompose 
the Green tensor in the following way: 



G(r,r» 



G (r,r',u) + G R (r, r',w) r>R 
G T (f,f\ui) r<R 



(2) 



where, if the observation point r is also outside the cylinder, 
we have the sum of a direct (or vacuum) and a reflected contri- 
bution, whereas inside the cylinder we get only a transmitted 
term, r is the distance from the wire axis. To fulfill the bound- 
ary conditions, we expand the Green tensor in cylindrical har- 
monics |9J|- The cylindrical harmonic vector wave functions 
are given by 



I Mo, n (fez, f) J 



Vx[z„(fc ro 1 r)(cosii0)e' fczZ i] 
Vx[z n (k r(l 1 r)(sinn0)e ifc z z z] 

r^— VxM c „(k*,r) 
■jJ— VxA? „<k z ,r) 



(3) 



Finally, we have for the transmitted part 

G T (r,r» = 

^ / dt L — js— 



(6) 



[^ T Me„(^, r) + B T No n {k z ,r)) M$(-k z ,?') 

+ (c T Ne n (k z ,r) + D T Mo n (k z , f)) N^(-k z , r% 

The reflection and transmission coefficients 
Ar.t, Br,t,Cr,t, Dr^t can be specified by imposing 
the boundary conditions at the surface of the cylinder 

(r = R): 

f x (5(r, f )r=fl- = f x G(r,r') r= _R+ (7) 
fxVxg(f,f'U- = fx V xG(r,r') r=R+ . (8) 

Solving the above equations, we get the full Green tensor. 
A. longitudinal modes 



with k = ~, ki = v/eTfco and kf. Q1 



fcg i — k\ where 



we have 0, 1 for outside and inside the cylinder and the sub- 
scripts e and o stand for even and odd, respectively. Note that 
k\ is complex due to the non-vanishing imaginary part of the 
electric permittivity e" = Im[e]. The radial part Z n (x) has 
to be either replaced by the Bessel function J n (x), or, when 
the "(1)" appears as a superscript on M or N, by Hn\x), 
i.e. the Hankel function of the first kind. Note, further- 
more, that the tensorial products between the even and odd 
cylindrical vector wave functions are defined according to 

Me n N^ = M e ,n/yi,n + M 0) „^vi,n and all other combina- 
tions similarly. For better readability, we have omitted the 
tensorial product symbol between them. 

This gives for the direct term of the Green tensor 



G Q (f,f',uj) = 
rrS(f— f') i 



dk z ^ 

-oo n=0 'To 



[M. n (k„T f )M$(-k z ,?') + N %n (k z ,r)N^(-k z ,f) 

where the upper (lower) line holds for r > r' (r < r'). For 
the reflected part we find 

G R {r,r',uj) = 



(1), 



— / dk z y 2 —^ 



[(A R M^{k z ,r) + B R N^(k z ,r)) M$(-k z ,r') 
+ (c R N^(k z ,r) + D R M$(k z ,r)) N§>{-k„ r% 



Let us first discuss the solution of the boundary equa- 
tions in the absence of material losses, i.e. for e" = 0. 
Writing up (0 and ©, one finds a system of linear alge- 
braic equations for the transmission and reflection coefficients 
A r ,t, B Rt T, C r .t, D R ,t- One can rewrite these as a vector, 
containing the coefficients, multiplied by a matrix. In order 
for the equations to have a nontrivial solution infinitely far 
from the source (i.e where we can neglect the contribution 
from Go), the determinant of the matrix has to vanish. From 
this condition we obtain a mode equation for each harmonic 
order n. These mode equations are equivalent to those ob- 
tained by different methods in 0]] : 



R 2 



1 



1 

K r 



1 J^(k ri R) 1 H'(k ra R) 



(9) 



(4) 



My n > (k z ,7)Me n (-k z ,r') + Ny n > {k z ,f)N %n (-k z , f) 



(5) 



k ri J n {k ri R) k ro H n (k ro R^) 
fcf J^K.R) k 2 H^(k ro R) 

fori Jn{kriR) k rQ B n {k rQ R^ 

Fixing the frequency oj, we get an equation for the longitudi- 
nal component of the wavevector parallel to the wire, the roots 
giving the allowed values k z = k z (uj; n) for each transversal 
mode n. These are the surface plasmons. In FigQ]the allowed 
values of k z are plotted as function of wire radius R for the 
different cylindrical harmonics. One recognizes that with the 
shrinking of R, all solutions belonging to n ^ modes have 
a cutoff and only the one belonging to n = persists. The 
corresponding value of k z diverges as R — > . 

When calculating the reflected and transmitted Green ten- 
sors, as seen in © and (0, the terms inside the integral and 
the sum on their right hand sides have a singularity exactly at 
k z = k z (uj;n). This means that there is a contribution to the 
response of the nano-wire system that originates from the n-th 
plasmonic mode. 
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FIG. 1: k z components of the guided modes of a nanowire (e = 
—75), as the function of the wire radius 7? and cylindrical harmonic 
order n. When the radius decreases, all modes with a non-zero n 
have a cutoff, whereas the n = mode has an increasing k z , that 
is, an icreasingly strong confinement. R is scaled with the vacuum 
radiation wavelength Aq. 
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f.TQ(n o) (ijt A ^ ta,u\ fc z )fj , in case of a single-mode, 
75 + 0.6i). f is the unit vector in the radial direc- 



FIG. 2: Im 

lossy wire (e 

tion. The broadened plasmonic resonance peaks show up at around 
the plasmonic longitudinal wavenumbers ±fc p ; (predicted by ((9]l), in 
the evanescent region (k p i > ko). The distances and wave numbers 
are scaled with the vacuum radiation wavelength Aq. 



B. width of the plasmon resonances 

Taking into account the medium absorption, i.e. a small but 
finite value of e", the discrete cylindrical modes turn into res- 
onances. This means that in case of losses there is no longer 
a single, well-defined k z but continuously many k z values 
peaked around k z (uj; n) that contribute to the plasmon mode. 
The consequences can be drawn much clearer if the wire ra- 
dius is much smaller than the optical wavelength, because then 
we can interpret the propagation of the plasmon mode as a ID 
problem, looking at a single mode only. The material losses 
induce a distribution of components with different k z values in 
the propagating plasmonic mode, each acquiring a phase fac- 
tor exp(ik z z) while propagating a distance z. This, in turn, 
results in a dephasing of the components upon propagation, 
which describes the losses. 

Fig. [2] shows the shape of the resonance belonging to the 
n = cylindrical mode around k z = ±k z (uj;n), the value 
corresponding to the lossless case. The lineshape is well ap- 
proximated with a Lorentzian function. The HWHM of the 
resonance peak scales linearly with e", as one would expect. 
In Fig|3]the value of the HWHM is depicted, which assumes 
larger and larger values as the wire radius R gets smaller. One 
recognizes a transition from a slower R~ x l 2 to a faster i? -3 / 2 
dependence of the width of the plasmon resonances at around 

10-%. 

So, in case of e" ^ 0, decreasing the wire radius is accom- 
panied not only by a decreasing transverse mode area, but also 
by increasing propagation losses. The two effects are con- 
nected because the smaller the transverse mode area gets, the 
larger part of the propagating field will be concentrated inside 
the metal which induces growing propagation losses. We will 




FIG. 3: HWHM of the plasmonic resonance peaks shown in Fig(2] 
as a function of wire radius R, in a log-log plot. The width and the 
radius are scaled by vacuum radiation wavelength Ao. The electric 
permittivity is e — — 75 + 0.6i. As the wire gets thinner, the broaden- 
ing of the resonance increases, i.e., the propagation losses increases. 



see later that for a decreasing value of R the detrimental ef- 
fects of propagation losses will overcome the positive effect 
of an increased coupling, resulting in an optimal wire radius. 
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III. ATOM-PLASMON COUPLING 
A. general framework 

A two-level atom with an electric dipole transition of 
strength d couples to the guided plasmon waves via the quan- 
tized electromagnetic field E(f A ) at the position of the atom. 
The interaction Hamiltonian in the dipole and rotating wave- 
approximations in the Schrodinger picture reads 



(+) 



where 



Ht = E (r A ) +h.a 



E {+) {r) = / du)E{r,u). 
Jo 



(10) 



(ID 



is the positive frequency part of the electric field operator in 
the presence of the nano-wire and = \e)(g\ is the atomic 
raising operator between the ground state \g) and the excited 
state |e) of the atom. The negative frequency counterpart of 

E\+> (r) is defined as 



E { - ] {r) = / dw#(r». 



(12) 



Following the treatment of B10I1 . the electric field operator 



h w 2 



7re C 



E(?,u)=iJ ¥ / d 3 fVe"(r', W )G(r,f»/ u (f') 



(13) 

can be expressed in terms of elementary excitations f u (r) and 
/J, (r) fulfilling bosonic commutation relations 



of an atom with the different modes of the electromagnetic 
field is inversely proportional to the square root of the ef- 
fective mode volume. That also means that the probability 
of spontaneous emission into modes with the smallest effec- 
tive mode volume will be the largest. Additionally, slow 
group velocity of the given mode also enhances the interac- 
tion strength. 

In the following we will take a look at the spontaneous 
emission rate of a two-level atom near the surface of a 
nanowire. We have seen before that it is possible to have 
a wire with a single, strongly confined mode at a given fre- 
quency, that is, a mode with a very small effective cross sec- 
tion area. Thus, we expect that if we place the atom close 
enough to the wire, its spontaneous emission rate will increase 
considerably | To ascertain this, we write up the Heisen- 
berg equations of motion for the field and atomic operators. 

'a z = ja f E (+) (r A ,t)d + H.c. (17) 

*t = ilua^ + l -E ( -\f A ,t)da z (18) 

where the operator of the electric field is given by dot . The 
equation of motion of the electric field, respectively the ele- 
mentary excitations / reads 



fu>(r) = -iufw(r) + 

+ u^_ j e"(r,uj) 
c 2 y hne Q 



G*(r,r A ,uj)da (19) 



To eliminate the field, we formally integrate ( fT9b and substi- 
tute it into ( fTTb and ( fT~8b . Employing a Markov approximation 
we get rid of the time integral, and using the Green tensor 
property (Hi) 



k fl>. = <W~ r')S(w - J), (14) 



o. 



(15) 



#<f = 








H A = 


-hu> A a z 



Here, e" is the imaginary part of the electric permittivity. The 
interactionless Hamiltonians of the combined field-matter sys- 
tem and a single, two-level atom have the simple form 



(16) 



where uo A is the atomic transition frequency and a z = \e)(e\- 

\g)(g\. 



B. spontaneous emission and level shift of an atom near a 
nanowire 

It is well known that a dielectric body near an atom influ- 
ences its spontaneous emisson rate. The interaction strength 



= Im G(fi,f 2 ,u}) 
we obtain the effective equations for a z and 

o z = -r tot (i + a z 



21 "-tK(+) 



a^EY 1 d + H.c. 

h free 



(20) 



(21) 



i(oj A - <M - -Ttot 



\E^)Jo z (22) 



where E^ ee = E'^ > ee {f Al t) is the free part of the field pro- 
portional to the free part of / and /'. The total decay rate and 
the Lamb shift T to t and Su> can then be expressed in terms of 
the electromagnetic Green tensor 



Ttot — 
Slu = 



2uj A didj _ 
fQ c 2 Im [ G vi r A,r A ,uj A )] 

did 



(23) 



I ^^^M, (24) 



OJ — 0J A 
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FIG. 4: Spontaneous decay rate of a single emitter relative to the vac- 
uum spontaneous decay rate as a function of emitter - wire axis dis- 
tance ta and wire radius R, both scaled by vacuum radiation wave- 
length Ao. We use e = —75 + 0.6i. Because of the strong interaction 
of the emitter with the wire eigenmodes, there is a considerable in- 
crease in the spontaneous decay rate at small radii and emitter-wire 
distances. Since Ttot — > oo for ta — > A, we have introduced a cutoff 
at a small ta — R for better visibility. 

Fig. E] shows the results of the numerical simulation for Ttot 
for different atom - wire axis distances and wire radii. Ac- 
cording to these results , the interaction of the atom with the 
surface plasmon mode gives rise to a significant Purcell ef- 
fect, enhancing the spontaneous emission rate by two orders 
of magnitude. Thus, the atom interacts with the surface plas- 
mon mode in the strong coupling regime. It is important to 
note that the interaction of the atom with the wire strongly de- 
pends on the direction of the polarization vector d. The cou- 
pling to the plasmons is strongest when d points in the radial 
direction. In the simulation we used this configuration. 

We have yet to justify the Markov approximation. In Fig. [5] 
we show Ttot, as a function of different atom-wire distances 
and frequency, taking a linear dispersion (this is a good ap- 
proximation at optical frequencies). The results show that Ttot 
has a very large width in frequency, much more than 0.5u>a, 
which means that the recurrence time (which is proportional to 
the inverse of the width) is several orders of magnitude smaller 
than the characteristic time scale of the atomic decay. This 
justifies the use of the Markov approximation 



C. sub-wavelength nano-wire 

Although with the Green tensor formalism one can handle 
the problem in an elegant way, it describes the interaction of 
one or many emitters with their surroundings in a compact 
form where contributions of different effects are difficult to 
distinguish. For example, when calculating the spontaneous 
emission of a single atom, we obtain the total emission rate 
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FIG. 5: Spontaneous decay rate of a single emitter relative to the 
vacuum spontaneous decay rate as a function of wave number k = 
ui/c and emitter - wire axis distance ta, in case of a lossy wire (e = 
— 75 + 0.6i) with radius R = O.OlAo. Both k and ta are scaled 
by vacuum radiation wavelength Ao. The large width in frequency 
verifies the applicability of the Markov approximation. Since Ttot — > 
oo for t a — ► R, we have introduced a cutoff at a small ta — R for 
better visibility. 



into all possible modes of the environment. Thus, on first 
glance, the result does not enable us to separately determine 
the rate of spontaneous emission into electromagnetic modes 
associated with the nano-wire plasmons and those into the ra- 
diative modes. In the following we will show, however, that 
with some further calculations we can still resolve the contri- 
bution of the different modes. We can define the following 
Hamiltonian: 

&! = Hf + HT (25) 

where we have separated the interaction into the evanescent 
modes (\k z \ > ko) and the traveling-wave part (\k z \ < ko) of 
the electromagnetic field. We do this because, as one can see 
in FigQ] the plasmons are evanescent in the radial direction 
(\k z \ > ko), never appearing in the traveling-wave part of the 
/c-spectrum. Moreover, for lossless media e" = 0, one can 
verify from the structure of lm[G^ n ' {ta, ta, k z )] that the 
surface plasmons are the only contribution in the evanescent 
region. In this way we can define traveling and evanescent cre- 
ation and annihilation operators a^ r,ev , a^'™ with the usual 
commutation relations, so that 

poo poo 

Hj = - du)g tr {u))ofctf - / dug ev (uj)aZ v ^ + h.a. 
Jo Jo 

(26) 

where g ev ' tr (uj) is the coupling strength to the evanescent and 
traveling modes, respectively, at frequency lo. is the atomic 
flip operator. Writing up the Heisenberg equation for a™' tr 
and a, and using that the traveling and evanescent operators 
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commute, we get 



(27) 



-iujA<j(t) 



--J du (<7 4 >)^ r + 9 eV {"YcC) *.(*)■ 

(28) 

Formally integrating d27T ) and substituting into d28l ) yields: 

a = — iu)Aa(£) 

oo 

cfco (^(^(O) + <7 e »a-(0)) e^S^f) 

+ ^ (l5 tr (^)| 2 + ls e »| 2 ) fdra{r)e^-^a z (t). 

(29) 



Note that in d29l there are no cross-terms because the trav- 
eling creation and annihilation operators commute with the 
evanescent ones. 

We can repeat the same procedure using the Green function 
formalism. For the equation of motion of a we obtain 



-ZU>AO'(t) — 



duji 




x d\'^(^)^G{r A ^,uj)f^)a z {t)OQ) 



By formal integration, we express the solution of (fT9b and sub- 
stitute it into (l3Qb . Using the Green function property (f20b for 
nonmagnetic materials we obtain the following expression: 

a = -iujAO - -&fr e JfA,t)da z {t) 
1 f x hw 2 



doj- 



rcFlm 



?i 2 Jo 7r£ oc 
x / dTa{ T )e-' lul(t - T) a z {t) 



G{r A ,r A ,u) 



(31) 



which is equivalent to d29l . Both the source and the observa- 
tion points are in vacuum so |fc| = fco = uj/c. Furthermore, 
we have seen in Section [TT] that we find the Green tensor by 
integrating an analytic, k z - and n-dependent expression over 
k z and summing it up over n. With this we can rewrite (l3~TT l 
which assumes the form 



<r = -iu A a - -E ( fr l e (r A ,t)da z (t) 

1 f°? H, ^ hio 2 
+—JdujJ dk z y j - 

ri 



\ 



/dTa{r)e 
-oo 



n=0 
— iuj(t — r) 



(32) 



Comparing (|29l and (1321 , separating the evanescent and 
traveling contributions, and exploiting the symmetry of 



Im 



r A ,r A ,to;, 



in k z we can identify 



ls e »l g 



poo °° 

dk z Y, 

J fen n 



hu 2 

ire c 2 



xd Im 



n=0 

\&^\r A ,r A ^k z ) 



d, (33) 



»| 5 



= 2 



/■fen oo 

/ dk *Yl 



hoj 2 

■KEqC 2 



xd Im 



(34) 



According to Fermi's golden rule, the square of the coupling 
strength to a given mode of the field is proportional to the 
spontaneous decay rate into that mode. Because of this and 
Eqs. ( 1331 and d34l l. we can define the decay rate into the plas- 
mon mode associated with the n-th cylindrical order: 



n=0 



dk, 



2u 2 
heoc 2 



(Flm 



(35) 



The above decomposition into traveling and plasmonic com- 
ponents is strictly valid only in the absence of losses. In case 
of a lossy wire, the plasmon resonances will not be singular- 
ities anymore but finite peaks with a finite width. Moreover, 
they will no longer be the only contribution in the evanes- 
cent region: apart from the resonance peaks, we observe a 
broad, low-amplitude background that is due to dissipative 
local circulating currents. As shown in Ref. ifTHl . e" ^ 
results in non-radiative losses which, in turn, have a contri- 
bution to \g tr (oj)\ 2 , and, because of the circulating currents, 
also to |g e "(w)| 2 . Our main interest lies in a thin, single- 
mode wire where we do not have higher-order plasmon modes 
(see FigQ]). This allows us to resolve the plasmon mode from 
the circulating surface currents to a good approximation. The 
evanescent part of the n > cylindrical orders contains only 
the circulating current background. As for the n = con- 
tribution, one can verify that there is no peak-sitting-on-a- 
background behaviour in Im [G < -™ =0 - ) (r A , r A , w; k z )], but in- 
stead we have a single peak, corresponding to the plasmon 
mode, with an almost perfect Lorentzian shape. So, in n = 
order we have no, or negligible circulating surface current 
contribution. Thus, in the single-mode case the decay rate 
into the plasmons can be calculated as 



T pl (uj) = 2 / dk 



2uj 2 
' he$c 2 



& n =°\r Al r A ^-k z ) 



d, 

(36) 



Plotting the emission rate into the plasmons normalized to 
the overall decay rate as a function of the atomic distance from 
the wire axis, as seen in Fig. [6] we observe that there is an 
optimum distance at which this value is maximal. This is in 
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FIG. 6: Spontaneous emission rate into the plasmons relative to the 
total spontaneous decay rate, as a function of source distance from 
the wire axis ta, scaled by the vacuum radiation wavelength Ao, in 
case of a wire of radius R — O.OlAo and e = — 75 + 0. 6i. There is an 
optimal ta at which the coupling of a single emitter to the plasmons 
is maximal. 



accordance with the result in Our model enables us to de- 
termine the origin of this effect. If the atom gets far from the 
wire, intuitively, the spontaneous emission rate into the plas- 
mons decreases rapidly, while the emission rate into the travel- 
ing modes stays more or less constant (practically zero reflec- 
tion from the wire and constant vacuum spontaneous emission 
rate) which results in the overall decrease of r p ;/r tot . How- 
ever, if the atom gets closer and closer to the wire, the circu- 
lating current contribution increases faster than the decay rate 
into the n = plasmon mode. This, again, leads to the overall 
decrease of the relative decay rate. 



IV. WIRE-MEDIATED INTERACTION OF EMITTERS 

A. sub- and super-radiant coupling mediated by the nano-wire 

If we place two emitters along the single-mode wire, both 
at the same distance from the surface, close enough to couple 
to the single surface plasmon mode, they will have a long- 
range interaction mediated by the plasmons. In the following, 
we will investigate this interaction. The rotating-wave Hamil- 
tonian of two identical, two-level atoms interacting with the 
environment in dipole approximation assumes the form 

H = J d s r J dujfkvfl(r)f u (f) + ^huj A {a zl + & z2 ) 
+ ( y -&\E {+) {f u t)d l ~&lE {+) {r 2 ,t)d 2 + H.c^ 

(37) 

where the indices 1 and 2 refer to the 1st and 2nd atom. Sim- 
ilarly to the treatment in fl2Tl . we can write up the master 
equation, trace out the reservoir and after applying Born and 
Markov approximations we get the equation of motion for the 



reduced, two-atom density matrix: 



El„/, + 
— \b\o k p + po\o k 



k,l=l 



(38) 



i 2J 8wk,i 
kd=i 



Sujkk are the single-atom Lamb shifts and 8u)i 2 is the radiative 
dipole-dipole shift. Similarly, is the single-atom decay 
rate T to t seen in Section lTlI Bl and Ti2 is a contribution which 
describes the effect of atoms upon each other. 



Tl2 



2u A d u d 2 



dud. 



2.1 . 



heoTT 



doj 



uj 2 Im [Gy(fi,r 2 ,w)] 



(39) 



(40) 



The one-atom flip operators have the form a k = \g) kk (e\, 
with \g)k being the ground state and \e) k the excited state 
of the fc-th atom. Since I^i = Ti 2 Eq. ( 1381 ) is diagonal in 
the basis of symmetrized and antisymmetrized states, namely, 
|ee>, \S) = (\ge) + \eg))/y/% \AS) = (\ge) - |e 5 ))/V2and 
\gg). Writing up the equations of motion for the populations 
of these states we get 



ree 

pil 

p% 



-WllPfe 

(r u + r 12 ) P £ 
(rn - v 12 )pii 



(r u + r 12 ) P | 
(r u - ri 2 )4| 



(r u + ri2)p| + (r u -ria)p[ 



AS 
AS- 



(41) 
(42) 
(43) 
(44) 



Additionally, from the equations of the coherences, we per- 
ceive that the resonance frequency of \AS) has been shifted 
downwards by &j 12 , whereas that of \S) has increased by 
the same amount. We used the following notation for the 
matrix elements of the two-particle density operator: pff = 
(gs\p\es). The magnitude of Ti2 will be responsible for the 
super- and subradiance of the symmetric and antisymmetric 
transitions. 

Fig. [7] and Fig. [8] show the value of the cross relaxation 
rate Ti 2 normalized to the single-atom decay rate Tn = T22, 
where we have varied the distance between the two atoms by 
several vacuum wavelengths while keeping the atom - wire 
distance constant, with metal losses switched off and on, re- 
spectively. We observe an oscillatory behaviour of Ti 2 be- 
tween ±Tii, where the period of oscillations matches exactly 
the longitudinal wavelength of the surface plasmon mode. In 
essence, this effect is caused by the interference of the radi- 
ation emitted by the two atoms. Depending on the distance, 
the two plasmon contributions emitted by the two atoms while 
being in the symmetric (or antisymmetric) state may either in- 
terfere constructively or destructively, giving rise to enhanced 
or suppressed decay rates, respectively. This means that vary- 
ing the inter-atomic distance, we can induce superradiance for 
the symmetric transitions while suppressing spontaneous de- 
cay for the antisymmetric ones, or vice versa. Besides, one 
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FIG. 7: Variation of Ti2 /Fn as a function of the distance d between 
two emitters, in a lossless, single-mode wire of radius R — O.OlAo. 
The emitter - wire axis distance is ta = 0.015Ao. The oscillations 
result in the alternation between super-and subradiance of the sym- 
metric and antisymmetric atomic transitions. All distances are scaled 
by the vacuum radiation wavelength Ao. 
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FIG. 8: Variation of Ti2 /Tn as a function of the distance d between 
two emitters, in a single-mode, lossy wire (e — — 75 + 0. 6i) of radius 
R — O.OlAo. The emitter - wire axis distance is ta = 0.015Ao- The 
oscillation is damped because of the dissipation of plasmons upon 
propagation. All distances are scaled by the vacuum radiation wave- 
length Ao . 

can see in Fig. |7]that the maximum value of Ti2 is very close 
to Tn which means an efficient enhancement/suppression of 
sponateous decay of the corresponding transitions. Therefore, 
the atoms are indeed strongly coupled to the guided mode - 
not only is the overall decay rate greatly enhanced but almost 
the whole of the radiated energy is carried by the plasmons. 
If there are losses in the metal, however, the oscillations are 
damped. 

Finite material losses, i.e. e" 7^ give rise to another effect 
related to the discussion in Section IIIICI If we put an emit- 
ter too close to the wire surface, non-radiative losses become 
increasingly important. Thus although the coupling strength 



to the plasmons increases, the non-radiative channel will take 
over. Therefore, as Fig. [9] shows, there is, again, an optimal 
distance of the emitters, when the contrast of the super- and 
subradiance is maximal. However, the optimum distance is 
slightly different from the one found in Section MI CI Since 
we have two, strongly interacting emitters, the decay rate of 
the two-atom system into the plasmon mode will change dif- 
ferently as compared to the single-atom decay rate. Because 
of this, one can observe that the maximum in Fig.|9]is slightly 
shifted compared to that of Fig. [6] 

The efficiency of the coupling between two emitters 
separated by a fixed distance d will depend on several factors. 
In order to maximize the single-atom plasmon coupling 
the atom should be placed as close as possible to the wire 
surface. Since at some point non-radiative losses of the single 
atom increase more rapidly when approaching the surface 
there is an optimum atom-wire distance for a given wire 
radius R. On the other hand the propagating plasmons also 
experience absorption originating from the same material 
losses. These absorption losses depend on the structure of the 
mode function and thus also on the wire radius R. Wehave 
examined whether for a fixed distance d along the wire there 
is a wire radius allowing an optimal compromise between 
these two effects. In Fig. [10] we plot the enveloping curve 
of the maxima of Fi^/Tii as a function of inter-atomic 
distance, in case of three different wire radii. In each case 
we used a wire-emitter distance that allowed an optimal 
atom-plasmon coupling. One can see that at a smaller 
radius R, the achievable maximum of T^/Tn is higher. 
However, as the propagation distance increases, the value of 
the envelope decreases faster compared to the case of a larger 
R. According to this, if the emitters are close to each other, 
it is better to have a wire with a smaller radius. On the other 
hand, if they are further apart the optimal radius would be a 
larger one. There is thus an optimal wire radius for a given 
inter-atomic distance, but there is no single optimal R for all 
cases. 

As a possible application, in the following we are going to 
analyze the realization of a deterministic quantum phase gate 
between two lambda atoms coupled by a single surface plas- 
mon mode of the nanowire, as shown in Fig QT| We assume 
that only the \g) — > \e) transition of each emitter is coupled to 
the single guided mode of the wire. If the two emitters are sev- 
eral wavelengths apart they can individually be addressed by 
lasers coupling the \g) — > |e) transition. The phase gate works 
by exploiting the large difference in the decay rates between 
the super- and sub-radian states shown in Fig. [8] Assume first 
that we can ignore the decay of the excited level |e) in the 
I. 9) ~ \ e ) two level system. In this case it is well known that a 
2n pulse, forcing the system to do a full Rabi oscillation, gives 
an additional it phase to the atomic system. On the other hand 
if the decay rate of the excited state r eg is much stronger than 
the resonant Rabi frequency O <C T, the driving field can- 
not induce Rabi oscillations between the two levels, and the 
drive merely introduces a scattering rate oc il 2 /T eg , which 
vanishes in the limit of weak driving or strong decay. The 
idea behind the phase gate is that due to the subradiant effect 
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FIG. 9: Magnitude of an extremal value of F^/Tn as the function 
of the atom - wire axis distance ta, scaled by the vacuum radiation 
wavelength Ao, at wire radius R = O.OlAo and inter-atomic distance 
d — l.OlAo. We use e = — 75 + 0.6i. There is an optimal ta where 
the coupling to the guided mode is maximal, resulting in a minimal 
damping of the oscillations of Ti2 /I\i . 
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FIG. 1 1 : Realization of a deterministic quantum phase gate, by cou- 
pling the |e) — \g) of two lambda atoms to a single surface plasmon 
mode, inducing superradiance in the two-atom system, and applying 
external classical 2n pulses. 
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FIG. 12: Effective level scheme of two, externally driven, interacting 
3-level lambda atoms, seen in Fig. QT] fis and Has are, respec- 
tively, the symmetric and antisymmetric superpositions of fii and 
0,2- Because of the strong coupling to a single plasmonic mode, a 
high-contrast superradiance effect is present in the system, inducing 
a slow decay in the symmetric transitions and a fast one in the anti- 
symmetric transition. 



FIG. 10: Maximal values of Ti^/Tii, as a function of propaga- 
tion distance d, for three different single-mode wire radii. We use 
e" — — 75 + 0.6i. The lines are exponential fits to the T12 /Tn peak 
values, represented by the markers. As R decreases, the propaga- 
tion losses increase but the largest achievable (ri2/rn) m aa; value 
also gets higher. Thus, at a given inter-emitter distance there exists a 
single R that ensures a maximal strength of interaction between the 
emitters along the wire. 



we saw above, the interaction of a nearby emitter can change 
the excited state from decaying into non-decaying thereby al- 
lowing the first emitter to pick up a it phase shift conditioned 
on the other other . 

Specifically, consider the level scheme of the two-atom 
system shown in Fig [12] where we use the symmetric- 
antisymmetric superposition for the transitions coupled to the 
wire. The external, classical fields with Rabi frequencies Oi 
and give rise to coupling strengths to the symmetric or an- 



tisymmetric states given by 



r»s = -^=(fii + n 2 ) 

n AS = - n 2 ) 



(45) 
(46) 



Assume now that the distance between the atoms is such 
that the symmetric state \S) is subradiant, whereas the anti- 
symmetric state \AS) is superradiant Fas ^ Ts and that 
= ^1 corresponding to Has = and Q$ = \/2f2i. 
Choosing the driving strengths £ls to be in between the two 
decay rates 



T A s > Ms > r s 



(47) 



ensures that transition to the double excited state |ee) is 
blocked by the strong decay Fas, whereas we can perform a 
2ir Rabi oscillation on the two level system given by \gg) and 
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\S) (see Fig. [T2| ). This thus ensures that we achieve a phase 
change of it on the \gg) state. Starting from the states \gs) and 
\sg) there is essentially only a single atom interaction with the 
wire. The excited state therefore always have a fast decay such 
that the driving is too weak to excite the atoms. Furthermore 
the state \ss) is completely unaffected by the classical light 
pulses and in the ideal limit the entire process will therefore 
have the truth table 

\ss) — > \ss) 
\sg) -> \sg) 

\gs) -> \gs) 

I.9.9) -\gg). (48) 

In reality the finite ratio r^s/rs will limit how well one 
can fulfill the condition in Eq. d47l i and this will limit the fi- 
delity F of the gate. Choosing the Rabi frequency too small 
will lead a decay from the excited state during the 2tt pulse 
on the \gg) — \S) transition. This will result in an error with 
a probability cx rsT pu i se oc r^r/fii, where the pulse dura- 
tion Tp U i so is set by T pu i so ~ 1/fii. Choosing the driving too 
strong will allow scattering from the transitions which are sup- 
posed to be blocked by the fast decay of the excited states with 
decay rates Ts oc T eg . This introduces an error probability 
oc T pu i se f2f /r eg oc Hi /r eg resulting in a total imperfection 

l-F^^L + ll. (49) 

1 eg "1 

Minimizing d49l with respect to fii, we find that the minimal 
imperfection scales as 

l-fopt~A&. (50) 

y i eg 

To confirm this rough scaling analysis we show in Fig [13] 
the result of a full numerical simulation of the density matrix 
equation for the system. In the figure we show the fidelity F = 
(V'ideallHV'ideal) optimized over the driving strength fii as a 
function of the ratio of the decay rates rs/r eg . Here IV'ideai) 
is a maximally entangled state created using the phase gate 
combined with ideal ir/2 pulses on the \g) — \s) transitions. 
The dashed line shows the scaling derived above 1 — F = 
c-v/T s /T eg , with the constant c chosen to match the behavior 
as F -+ 1. 

The analysis above has not been been tied to any particular 
system, but applies to to any system where there is a large 
difference between the super- and sub-radiant states. We now 
turn to the specific implementation in terms of atoms coupled 
through the plasmon. As we saw above, the subradiant state 
do not have a completely vanishing decay rate and accordingly 
this will limit the fidelity of the gate. In Fig. [14] we show the 
optimal fidelity of the gate operation, optimized by varying 
f2s, as a function of the atom - wire axis distance ta- For the 
simulation, we used R = 0.003Ao and d = 0.08Ao 

The results show a maximum fidelity of 80% for the param- 
eters above. There are two major detrimental effects which 
reduce the superradiance-subradiance contrast and, conse- 
quently, the gate fidelity, namely, the coupling to free space 




ryr 

S eg 

FIG. 13: Fidelity of a maximally entangled state created by the phase 
gate optimized with respect to the drive strength. The phase gate 
fidelity error scales roughly as y/Ts/T eg (dashed line) in the limit 
F-tl. 
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FIG. 14: Optimal fidelity of a maximally entangled state generated 
with the phase gate as a function of the atom - wire axis distance 
ta- ta is scaled with the vacuum radiation wavelength Ao. The two- 
atom states \ss), \sg), \gs) and \gg) are initially equally populated. 
We used parameters R = 0.003A , d = O.O8A0, e = -75 + 0.6i in 
the simulation. 

and the wire losses. Due to the free space coupling, a small 
amount of the radiated energy will be scattered into the far- 
field region, resulting in a non-zero contribution to Tg. The 
wire losses induce traveling-wave and circulating current dis- 
sipation - leading to an additional, local loss of energy. This 
also contributes to Ts- It is possible to suppress the sponta- 
neous emission into free space, for example, by placing the 
atoms into a photonic bandgap material. This would consid- 
erably enhance the gate fidelity only if the free-space contri- 
bution to Ts is dominant with respect to the wire-loss con- 
tribution. Fig. [15] shows the ratio Yg lre /Tg ree as a function 
of de" /Ao, for distances d corresponding to the minima of 
Ti2/rn. As one can see, to increase the inter-emitter dis- 
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FIG. 15: Ratio of wire-loss-induced and free-space-loss-induced 
contributions of Fs, as a function of e"d/Xo where we restrict d to 
take up values of the distances belonging to the minima of ri2/rn 
(see FigH. In the region of r^ ire /r f / ee < 1 the gate fidelity could 
be considerably increased by suppressing the spontaneous radiation 
into free space. 

tance d while keeping YJ ire /T f s ree < 1, we need to decrease 
e". For the given parameters, this means that, for instance, 
for d ps 2Ao we have to have e" ps 0.1. Another means of 
increasing the fidelity could be to build a cavity for the plas- 
mon, e.g., by using a ring shaped wire, but this will depend on 
the losses of such a cavity mode, e.g., due to radiation. A full 
investigation of this is beyond the scope of this paper. 



B. realization of spin boson model 

The plasmon mediated interaction between spins along a 
nano-wire could also be used to implement interesting and im- 
portant many -body models such as the spin boson model lfl3ll . 
in which a chain of spins is coupled to a one-dimensional con- 
tinuum of bosonic modes through an interaction of the type 

H SB = jf du £ (a, + aj) (g(uj)a u + g*(uj)at) . (51) 

Different from the atom-plasmon coupling ( fTOb discussed 
throughout this paper, this Hamiltonian contains counter ro- 
tating terms of the form a' a* and ad. Typically these terms 
can be neglected for electromagnetic modes oscillating at op- 
tical frequencies. However it is possible to compensate the 
fast oscillations by using effective two-photon Raman tran- 
sitions rather than single-photon ones [14]. Consider e.g. 
the F = 1/2- F = 1/2 coupling scheme shown in FigfT6l 
where the lower state \g) is coupled to another state \s) in 
the ground-state manifold via a two-photon Raman transition 
through an excited state |e+) involving a right cicular polar- 
ized external drive field of Rabi-frequency Q + and the linear 
polarized quantized plasmon field E. In addition state | s) is 
coupled through another excited state |e_) via a two-photon 




FIG. 16: Raman coupling scheme for the realization of a spin bo- 
son model. In the rotating wave approximationfhe atom can make a 
transition from \g) to \s) either by emitting a photon into E via the 
transition \g) — > |e+) — ¥ \s) or by absorbing a photon from E via 
the transition \g) — y |e_) — s> \s). 

Raman transition involving a left circular polarized external 
drive field il„ and again the linear polarized plasmon field 
E. It is assumed that the Raman transitions are in two-photon 
resonance. On the other hand the single-photon transitions 
shall have a large detuning A± such that the excited states 
can be adiabatically eliminated. For each single-photon tran- 
sition the rotating-wave approximation is well justified since 
we are considering optical fields. 

Hj = J2{ m *+*i+, g + n -ti-,s + h.a.} (52) 

3 

j 

Here M „ = are the atomic flip operators between 

states \v) and \fi) of the jth atom at position Zj, and d± are 

the vector dipole moments of the transitions |s) |e+) and 

\g) |e_) respectively. Adiabatically eliminating the ex- 
cited states then leads besides some ac-Stark shifts to the spin- 
boson interaction 

Hon - ^i^+^-^J^ (z 3 )+h.a) . 

j + 

(53) 

V. SUMMARY 

In the present paper we have analyzed the coupling of 
quantum dipole oscillators to the plasmon modes of an infi- 
nite cylindrical nano-wire by means of a Green function ap- 
proach taking into account wire losses. Placing a single emit- 
ter close to the surface of a cylindrical metallic wire with 
sub-wavelength radius, there is a strong interaction between 
the emitter and the traveling surface plasmon modes of the 
metal. The magnitude of the interaction becomes apparent in 
a substantial increase of the total spontaneous decay rate go- 
ing up to several hundred times the vacuum value. Decreas- 
ing the wire radius, only the basic cylindrical mode survives 
with a rapidly increasing radial confinement, allowing for a 
single-mode wire and an increasingly strong emitter-plasmon 
interaction. We have calculated the spectrum of the plasmon 
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modes and analyzed their broadening in the presence of metal 
losses (e" 7^ 0) for different wire radii R. For fixed e", the 
broadening and so the propagation losses increase with de- 
creasing wire radius since a growing part of the field is con- 
centrated inside the wire. For a sub-wavelength nano-wire the 
dependence of the resonance width on the wire radius turns 
from an i?^ 1 / 2 to an i?~ 3 / 2 behaviour. We derived explicit 
expressions for the coupling strength of an atom to the plas- 
mon modes and discussed its relation to non-plasmonic cou- 
plings in dependence on the wire parameter and the atom-wire 
distance. For non-vanishing material losses there is an opti- 
mum emitter-wire distance for which the decay rate into the 
plasmons relative to the total decay rate is maximal. We also 
discussed the plasmon-mediated coupling between two quan- 
tum dipoles. Putting two emitters close to the wire results 
in Dicke-sub- and superradiance even when the separation of 
the emitters along the wire exceeds the vacuum wavelength 
by an order of magnitude. Depending on the precise distance 
between the emitters, the decay rate of the symmetric or the 
anti-symmetric Dicke state is enhanced or suppressed. Simi- 
lar to the one-atom case, due to the finite material losses there 
is an optimum emitter-surface distance yielding the maximal 
coupling between the atoms. Also, for a given inter-atomic 
distance there exists an optimal wire radius for which the cou- 
pling between the emitters is maximal. The possibility to con- 
trol the symmetric and anti-symmetric decay rates allows to 
build two-qubit quantum gates using the setup. We proposed 
and analyzed a scheme for a deterministic two-qubit phase 
gate. The achievable gate fidelity for realistic parameters is 



on the order of 80%. This can be attributed to two loss mech- 
anisms: wire absorption of plasmons and spontaneous emis- 
sion into free-space radiation modes. For small values of e" 
the fidelity could be improved by embedding the atoms in a 
photonic band-gap material, suppressing the free-space spon- 
taneous emission. Furthermore, introducing a wire ring in- 
stead of an infinite cylinder, the effect from the loss of ex- 
citation along the wire may be reduced while retaining the 
strong coupling. Additionally, a geometry-independent anal- 
ysis shows that, asymptotically, the gate fidelity error scales as 
^jTs/T eg , i.e., the second root of the ratio between the sym- 
metric (subradiant) and one-atom decay rate, as F — > 1. Fi- 
nally atoms coupled to plasmonic nano-wires may also be an 
interesting system for the realization of interesting spin mod- 
els such as the spin-boson model. 



VI. ACKNOWLEDGEMENTS 

David Dzsotjan acknowledges financial support by the 
EMALI Marie-Curie Network and by the Research Fund of 
the Hungarian Academy of Sciences (OTKA) under contract 
No. 78112. Anders S. S0rensen acknowledges the support of 
the Villum Kann Rasmussen foundation and the Danish Na- 
tional Science Research Foundation. The authors would like 
to thank R Rabl for discussions and his valuable contribution 
to Section III B. 



[1] D. E. Chang, A. S. S0rensen, P. R. Hemmer, M. D. Lukin, Phys. 

Rev. Lett. 97, 053002 (2006); D. E. Chang, A. S. S0rensen, P. 

R. Hemmer, M. D. Lukin, Phys. Rev. B 76, 035420 (2007) 
[2] J. Takahara et al, Opt. Lett. 22, 475-477 (1997) 
[3] V. V. Klimov, M. Ducloy, V. S. Letokhov, Chem. Phys. Lett. 

358, 192-198 (2002); 1. 1. Smolyaninov, J. Elliott, A. V. Zayats, 

C. C. Davis, Phys. Rev. Lett. 94, 057401 (2005) 
[4] K. Kneipp et al, Phys. Rev. Lett. 78, 1667-1670 (1997); S. Nie 

and S. R. Emory, Science 275, 1102-1106 (1997); S. Kuhn, U. 

Hakanson, L. Rogobete, V. Sandoghdar, Phys. Rev. Lett. 97, 

017402 (2006) 

[5] D. E. Chang, A. S. S0rensen, E. A. Demler, M. D. Lukin, Na- 
ture Phys. 3, 807-812 (2007) 
[6] A. V. Akimov et al, Nature 450, 402-406 (2007) 
[7] A. L. Falk et al, Nature Phys. 5, 475-479 (2009) 



[8] R. Kolesov et al, Nature Phys. 5, 470-474 (2009) 
[9] L. W. Li et al, J. El. Magn. Waves Appl. 14, 961-985 (2000) 
[10] H. T. Dung et al, Phys. Rev. A 68, 043816 (2003); L. Knoll, 

S. Scheel, D.-G. Welsch, QED in Dispersing and Absorbing 

Media, In: Coherence and Statistics of Photons and Atoms, ed. 

by J. Pefina, Wiley (New York ,2001) 
[11] M. Fleischhauer, Phys. Rev. A 60, 2534 (1999) 
[12] J. Kastel and M. Fleischhauer, Phys. Rev. A 71, 011804(R) 

(2005); J. Kastel and M. Fleischhauer, Laser Phys. 15, 1, 135- 

145 (2005) 

[13] I. Bloch, J. Dalibard, W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008) 

[14] F. Dimer, B. Estienne, A. S. Parkins, H. J. Carmichael, Phys. 
Rev. A 75, 013804 (2007) 



